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Abstract 

Tree-level unitarity constraints on the masses of the Higgs bosons in the general 
Two Higgs Doublet Model (THDM) are studied. We first consider the case where 
the Higgs potential is invariant under a discrete symmetry transformation, and derive 
strong constraints on the mass of the lightest CP-even Higgs boson (Mh) as a function 
of tan (3. We then show that the inclusion of the discrete symmetry breaking term 
weakens the mass bounds considerably It is suggested that a measurement of Mh 
and tan (5 may enable discrimination between the two Higgs potentials. 
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1 Introduction 



The Minimal Standard Model (MSM) of electroweak interactions |IJ is in complete agree- 
ment with all precision experimental data (LEPII,Tevatron, SLD), with the notable excep- 
tion of neutrino oscillation experiments 0. In this minimal version there is one complex 
SU (2) i (8) £7(1) doublet which provides mass for the fermions and gauge bosons. After elec- 
troweak symmetry breaking three of the four real degrees of freedom initially present in the 
Higgs doublet become the longitudinal components of the gauge bosons W ± , Z, leaving 
one degree of freedom which manifests itself as a physical particle (<f)°) ||. The Higgs 
mass (Mtjfi) is not fixed by the model, although constraints on M^o can be obtained by 
making additional theoretical assumptions. These include the unitarity bound (M^o < 870 
GeV) [||] and the triviality bound || . So far no experimental information on the nature of 
the Higgs particle has been found, and from the negative searches in the Higgsstrahlung 
channel at LEPII the lower bound M^o > 107.7 GeV |6j has been derived. 

In recent years there has been growing interest in the study of extended Higgs sectors 
with more than one doublet [[/J. The simplest extension of the MSM is the Two Higgs 
Doublet Model (THDM), which is formed by adding an extra complex SU(2)l <g> U(1)y 
scalar doublet to the MSM Lagrangian. Motivations for such a structure include CP- 
violation in the Higgs sector, supersymmetry, and a possible solution to the cosmological 
domain wall problem ||. In particular, the Minimal Supersymmetric Standard Model 
(MSSM) takes the form of a constrained THDM. 

The most general THDM scalar potential which is renormalizable, gauge invariant 
and CP invariant depends on ten parameters, but such a potential can still break CP 
spontaneously |J. In order to ensure that tree-level flavour changing neutral currents are 
eliminated, a discrete symmetry ($« — > — $j, where is a scalar doublet) may be imposed 
on the lagrangian [[TIJ], which reduces the number of free parameters to 6. The resulting 
potential was considered in [IT|, and is referred to as Va in [O. We shall be concerned 



with the potential described in which is equivalent to Va plus a term which breaks the 
discrete symmetry (parametrized by A5) and contains 7 free parameters. Such a potential 
does not break CP spontaneously or explicitly |TT],(T3| provided that all the parameters 
are real. 

We note here that tree-level unitarity constraints for the THDM scalar potential Va 
were studied in [|14|, [15]. When deriving constraints from unitarity, |l4j considered only 



seven elastic scattering processes Si 62 — > Si £2 (where Si is a Higgs scalar) while |T5 



considered a larger scattering (S) matrix. In [15] , upper bounds on the Higgs masses were 
derived, in particular Mh < 410 GeV for tan/3 = 1, with the bound becoming stronger 
as tan {3 increases. We improve those studies for Va by including the full scalar S matrix 
which includes channels which were absent in [ 15 1 , and we also show graphically the strong 
correlation between and tan f3. 

To our knowledge such unitarity constraints have not been considered for the case of 
A5 7^ and this is the principal aim of this note. We shall see that the presence of a 



no n- zero A5 can significantly weaken the unitarity bounds found in [14, 15 



The paper is organized as follows. In Section 2 we give a short review of the THDM 
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potential and explain the unitarity approach we will be using. In Section 3 we present our 
numerical results for the cases of A5 = and 7^ 0, while Section 4 contains our conclusions. 



2. Scalar potential and unitarity constraint 



It has been shown Jl6j that the most general THDM scalar potential which is invariant 
under SU{2) L ® U(l)y and conserves CP is given by: 

V(*i, $2) = - vlf + A 2 (|$ 2 | 2 - vlf + A 3 ((|$ 1 | 2 - v\) + (|$ 2 | 2 - vl)f + 

A 4 (|$i| 2 |$ 2 | 2 - |$2"$ 2 | 2 ) + A 5 (i?e($+$ 2 ) - v lV2 f + A 6 [/m($+$ 2 )] 2 (2.1) 

where $1 and $ 2 have weak hypercharge Y=l, v% and f 2 are respectively the vacuum 
expectation values of $1 and $ 2 and the Aj's are real-valued parameters. 



Vi + 



V2 I 

This potential violates the discrete symmetry — > — softly by the dimension 2 term 
Asi?e($^$ 2 ) and has the same general structure of the scalar potential of the MSSM. One 
can prove easily that for A5 = the exact symmetry — > — $j is recovered. 

After electroweak symmetry breaking, the W and Z gauge bosons acquire masses given 
by m w = ^9 2y2 an d Trig = \{g 2 + g' 2 )v 2 , where g and g' are the SU(2) wea k and U(1)y gauge 
couplings and v 2 = v \ + v \. The combination v\ + v\ is thus fixed by the electroweak scale 
through v\ + v\ = (2\^2Gf)~ 1 , and we are left with 7 free parameters in eq. (P-l|), namely 
the (Ai)i=i ) ... ) 6's and tan/3 = Vijv\. Meanwhile, three of the eight degrees of freedom of the 
two Higgs doublets correspond to the 3 Goldstone bosons [G , G°) and the remaining five 
become physical Higgs bosons: H°, h° (CP-even), A (CP-odd) and H . Their masses 
are obtained as usual by the shift $j — > $j + [[F]. After generating the scalar masses in 
term of scalar parameters Aj, using straightforward algebra one can express all the A, as 
functions of the physical masses: 



\ y 2 \ y 2 \ y & ac a/ 2 2N ^5 

X4= 2^ mH± ' Ae = 2^ ' A3= 8^^ (m -- mJ " T (2 ' 2) 
Ai = [4ml + s 2 m 2 - - m 2 )] - ^(-1 + tan 2 0) (2.3) 



A 2 = a 2 2 [ s a m i + c a m h - s Q c Q tan/3(m| - m£)] - ^(-1 + — (2.4) 
8s^77% 4 tan z /3 

The angle (3 diagonalizes both the CP-odd and charged scalar mass matrices, leading to 
the physical states H and A . The angle a diagonalizes the CP-even mass matrix leading 
to the physical states H°, h°. 

We are free to take as 7 independent parameters (Xi)i=i t ...fi and tan/3 or equivalently 
the four scalar masses, tan/3, a and one of the A*. In what follows we will take A5 as a free 
parameter. 
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To constrain the scalar potential parameters one can demand that tree-level unitarity 
is preserved in a variety of scattering processes. This corresponds to the requirement that 
the J = partial waves (ao) for scalar-scalar and gauge boson scalar scattering satisfy 
| do | < 1/2 in the high-energy limit. At very high-energy, the equivalence theorem |T7[ 



states that the amplitude of a scattering process involving longitudinal gauge bosons 
may be approximated by the scalar amplitude in which gauge bosons are replaced by 
their corresponding Goldstone bosons G ±,a . We conclude that unitarity constraints can be 
implemented by solely considering pure scalar scattering. 

In very high energy collisions, it can be shown that the dominant contribution to the 
amplitude of the two-body scattering S1S2 — > S3S4 is the one which is mediated by the 
quartic coupling. Those contributions mediated by trilinear couplings are suppressed on 
dimensional grounds. Therefore the unitarity constraint \a \ < 1/2 reduces to the following 
constraint on the quartic coupling, \Q(SiS2S 3 S4)\ < 811. In what follows our attention will 
be focussed on the the quartic couplings. 

In order to derive the unitarity constraints on the scalar masses we will adopt the 



technique introduced in |y|. It has been shown in previous works |LJ| that the quartic 
scalar vertices written in terms of physical fields H ± , G ± , h°, H°, A and G°, are very 
complicated functions of A«, a and (3. However the quartic vertices (computed before 
electroweak symmetry breaking) written in terms of the non-physical fields (pf , hi and Zi 
(i=l,2) are considerably simpler expressions. The crucial point of [I5j is the fact that the 



S matrix expressed in terms of the physical fields (i.e. the mass eigenstate fields) can be 
transformed into an S matrix for the non-physical fields iff , hi and z% by making a unitarity 
transformation. The latter is relatively easy to compute from eq. [2.1| . Therefore the full set 
of scalar scattering processes can be expressed as an S matrix composed of 4 submatrices 
which do not couple with each other due to charge conservation and CP-invariance. The 
entries are the quartic couplings which mediate the scattering processes. 

The first submatrix corresponds to scatterings whose initial and final states are one of 
the following: ((pf^ ^Pii hx%2, h-iZx, Z1Z2, hih^)- Therefore one obtains a 6 x 6 matrix 
leading to the following 5 distinct eigenvalues: 

ei = 2A 3 - A 4 - y + ^A 6 

e 2 = 2A 3 + A 4 - y + ^A 6 

5 1 
/+ = 2A 3 — A 4 + -A5 — -A 6 

/- = 2A3 + A4 + — Ag — -A 6 

fx = h = 2A 3 + \h> + \\ (2.5) 

The second submatrix corresponds to scatterings with initial and final states one of the 
following: ((ftfi, ^2^2^ ^j 1 ? ^If)- Again we obtain a 6 x 6 matrix with the 
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following 6 eigenvalues: 



a± = 3(Ai + A 2 + 2A 3 ) ± t/9(Ai - A 2 ) 2 + (4A 3 + A 4 + ^(A 5 + A 6 )) 2 (2.6) 



b ± = Ax + A 2 + 2A 3 ± J(Xi - A 2 ) 2 + ~(-2A 4 + A 5 + A 6 ) 2 (2.7) 



c± = Ai + A 2 + 2A 3 ± y (Ai - A 2 ) 2 + ^(A 5 - A 6 ) 2 (2.8) 

The third submatrix corresponds to the basis: (hiZ\, /i 2 z 2 ). The 2x2 matrix possesses the 
eigenvalues d± and c±, with d± = c±. All the above eigenvalues agree with those found in 



T5[, up to a factor of 1/167T which we have factorised out. In our analysis we also include 



the two body scattering between the 8 charged states: h±ipi, h 2 f{, ziip{ , z 2 y?i , 



h2^>t > z i^P2 ) z 2¥>2 ■ Note that these channels were neglected in [15]. The 8x8 submatrix 
obtained from the above scattering processes contains many vanishing elements, and the 
8 eigenvalues are straightforward to obtain analytically. They read as follows: /_, e 2 , ft, 
c±, b± and pi, where 

Pl = 2(A 3 + A 4 ) - ^A 5 - l -\ (2.9) 

As one can see, these additional channels lead only to one extra eigenvalue, p\, although 
we shall see that this eigenvalue plays an important role in constraining M#± and Ma- 



3. Numerical results and discussion 

In this section we present our results for the unitarity constraints on the Higgs masses in 
the THDM. All the eigenvalues are constrained as follows: 

|o±l , \b±\ , \c±\ , |d±| , |/±| , |ei >2 | , |/i j2 | , bil < 8vr (3.1) 



In [15| the eigenvalues explicitly contain the factor 1/167T and so were constrained to be less 
than 1/2. In section 3.1 we consider the special case of A5 = 0, while section 3.2 presents 
results for the general case of A5 7^ 0. 

3.1 Case of As = 



For A5 = our potential is identical to those considered in [0, [1J|. We improve those 
analyses on two accounts: 

(i) We have considered extra scattering channels which leads to one more eigenvalue 
constraint, pi, as explained in Section 2. 

(ii) When finding the allowed parameter space of Higgs masses we simultaneously impose 



all the eigenvalue constraints. In |15j only the condition \a+\ <8ii was applied when 
deriving mass bounds. 
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M H ± 


M A 


M h 


M H 


AAN 


691 


695 


435 


638 


KKT 


860 


1220 


410 


700 



Table 1: A comparison of our bounds (AAN) and those of KKT [ Dj ] 
GeV and A 5 = 0. 



All masses are in 



In order to obtain the upper bounds on the Higgs masses allowed by the unitarity con- 
straints we vary all the Higgs masses and mixing angles randomly over a very large param- 
eter space. We confirm the result of [[TJJ which states that a + is comfortably the strongest 
individual eigenvalue constraint. However, the other eigenvalues impose important con- 
straints on Ma and M H ±. If only \a + \ < 8tt is imposed we can reproduce the upper bounds 
on the Higgs masses given in fi~5|| , in particular their main result of Mh < 410 GeV. When 
all eigenvalue bounds are applied simultaneously we find improved bounds on the Higgs 
masses, particularly for Ma and M#± . We note that the new eigenvalue constraint pi < 8tt 
(eq.2.9) plays a crucial role in determining the upper bound on Mjj±. We write in Table. 1 
our bounds (AAN) and those of P~5| , denoted by KKT. 

Note that the bounds given in Table. 1 are obtained for relatively small tan/3 (say 
tan/3 « 0.5). For large tan/? the bound is stronger, although for the case of A , H° and 
the tan j3 dependence is rather gentle. Of particular interest is the tan f3 dependence 
of the bound on M^ which will covered in the section 3.2. 



3.2 General case of A5 7^ 

We now consider A5 7^ which corresponds to the inclusion of the term which softly breaks 
the discrete symmetry. Such a term was neglected in the analyses of ]n, |T5|, and from 
perturbative constraints may take values | A5 1 < 8tt [19]. In the graphs which follow we 



do not impose the perturbative requirement |Aj| < 8tt for the remaining Aj. Imposing this 
condition only leads to minor changes in the numerical results which will be commented 
on when necessary. 

We plot in Fig.l the maximum value of Mh against tan/3 for increasing values of A5, 
imposing all the eigenvalue constraints simultaneously as done in Section 3.1. For the 
case of A5 = one finds a strong correlation, with larger tan/? requiring smaller Mh- For 
example, tan/3 > 7 corresponds to Mh < 100 GeV, which is the mass range already being 
probed by LEPII. However, h° in the THDM with Mh < 100 GeV is not guaranteed to 
be found at LEPII due to the suppression factor of sin 2 (/3 — a) for the main production 
process e + e _ — > h°Z. For the case of A5 = we find that values of tan j3 > 20 are strongly 
disfavoured since they easily violate one of the unitarity constraints. If A 5 7^ 0, Fig.l shows 
that for a given tan/3, the action of increasing A5 allows larger maximum values of Mh- 
For A5 = 15 one finds a horizontal line at Mh ~ 670 GeV, showing that the upper bound 
has been increased for all values of tan/3. In addition, large (> 30) values of tan/3 are 
allowed if A5 7^ 0, in contrast to the case of A5 = 0. For example, A5 = 1 comfortably 
permits values of tan f3 = 60. However, as pointed out in |2(J perturbative constraints on 
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Figure 1: Maximum in GeV as a function of tan/? for various values of A5. 



the Aj also restrict the allowed values of tan/3 in the THDM. Using the condition in 19 
which requires |Aj| < 8tt, we found that tan/3 > 30 is strongly disfavoured. In Fig. 2 we 
plot the maximum mass values (Ms) of all the Higgs bosons as a function of A5. Again 
all the Higgs masses and mixing angles have been varied randomly. We see that A5 = 
corresponds to the values in Table 1, while for A5 = 15 the upper bounds have been 
increased significantly. Including the perturbative requirement would lower the bounds on 
M A and M H ± by 10 -> 20 GeV. The tan/3 dependence of M s for S = H ± , A and H° is 
not very pronounced, and the maximum mass value may be obtained for both small and 
large tan (5. 

We note that the relaxation of the strong correlation between Mh and tan j3 with A5 7^ 
would in principle allow the possibility of distinguishing between the discrete symmetry 
conserving and violating potentials. If h° is discovered and the measured values of and 
tan j3 lie outside the rather constrained region for A5 = 0, this would signify A5 7^ and thus 
a soft breaking of the discrete symmetry. Possibilities of measuring tan /3 at high-energy 



e + e~ colliders have been considered in the context of the MSSM in [^T|]. Production and 
decay of H and A are particularly promising since the rates do not involve the mixing 
angle a. Much of the analysis of |21] is also valid in the THDM. 
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Figure 2: Maximum mass value of the Higgs scalars (Ms) in GeV as a function of A 5 . 

4. Conclusions 

We have derived upper limits on the masses of the Higgs bosons in the general Two Higgs 
Doublet Model (THDM) by requiring that unitarity is not violated in scalar scattering 
processes. We first considered the THDM scalar potential which is invariant under a 
discrete symmetry transformation and improved previous studies by including the complete 
set of scattering channels. Stronger constraints on the Higgs masses were derived. Of 
particular interest is the tan/3 dependence of the upper bound on M h , with larger tan/? 
requiring a lighter h° e.g. tan/3 > 7 implies Mh < 100 GeV. We then showed that the 
presence of the discrete symmetry breaking term parametrized by A5 may significantly 
weaken the upper bounds on the masses. In particular, the aforementioned correlation 
between tan (5 and maximum M h is relaxed. It was suggested that a measurement of tan (3 
and M h may allow discrimination between the two potentials. 
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